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In this paper, we develop a nonequilibrium theory for transient electron transport dynamics in 
nanostructures based on the Feynman- Vernon influence functional approach. We extend our previ- 
ous work on the exact master equation describing the non-Markovian electron dynamics in the double 
dot [Phys. Rev. B78, 235311 (2008)] to the nanostructures in which the energy levels of the central 
region, the couplings to the leads and the external biases applied to leads are all time-dependent. 
We then derive nonperturbatively the exact transient current in terms of the reduced density matrix 
within the same framework. This provides an exact non-linear response theory for quantum trans- 
port processes with back-reaction effect from the contacts, including the non-Markovian quantum 
relaxation and dephasing, being fully taken into account. The nonequilibrium steady-state transport 
theory based on the Schwinger-Keldysh nonequilibrium Green function technique can be recovered 
as a long time limit. For a simple application, we present the analytical and numerical results of 
transient dynamics for a resonant tunneling nanoscale device with a Lorentzian-type spectral density 
and ac bias voltages, where the non-Markovian memory structure and non-linear response to the 
time-dependent bias voltages in transport processes are demonstrated. 
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I. INTRODUCTION 

The investigation of quantum coherence dynamics far 
away from equilibrium in nanoscale electronic devices has 
been attracting much attentions in the past decade due 
to various applications in nanotechnology and quantum 
information processing. Experimentally it became possi- 
ble not only to directly manufacture structures but also 
to investigate their nonequilibrium quantum coherence 
properties under well controlled parameters. Theoret- 
ically, tremendous studies have been done focusing on the 
time-independent steady-state to analyze the nonequilib- 
rium transport properties. However, an increasing num- 
ber of work employing various theoretical techniques (see 
Refs. [8-24]) has recently occurred addressing the time- 
dependent electron dynamics far away from the equilib- 
rium. In fact, for practical applications, a full under- 
standing of nonequilibrium dynamics under external field 
controls, i.e., the time evolution of electrons in nano- 
electronics devices from some initial preparation toward 
any specifically designed state within an extremely short 
time, has been receiving more and more attentions. 

The prototypical nanostructure concerned in this pa- 
per is made by a gate-defined region on a semiconductor 
containing a few of discrete electronic states. Electrodes 
are implanted around this central region to control the 
electrons across it. Also gates are deposited to adjust 
the electronic states within the central area as well as 
their couplings to the surrounding electrodes. Most of 
the theoretical developments to nanoelectronics devices 
have been focusing on the understanding and prediction 



of transport dynamics in these devices. However, as a 
quantum device, in particular for quantum information 
processing,*^'' the big challenge is to understand and pre- 
dict not only how fast or slow a nanoelectronic device can 
turn on or off a current, but also how reliable and effi- 
cient the device can manipulate quantum coherence of 
the electron states through external bias and gate volt- 
age controls. This requires a nonequilibrium quantum 
device theory to analyze the time-dependent quantum 
transport intimately entangling with quantum coherence 
of electrons inside the device. It is the purpose of this 
paper to attempt to establish a nonequilibrium theory of 
transient dynamics for electron quantum states accompa- 
nied with electron quantum transport in a nanoelectronic 
device that can be feedbackly controlled through the non- 
linear response to the external time-dependent bias and 
gate voltage pulses. 

Physically the nanostructure we are going to study is 
typically an open quantum system in the sense that its 
central region has exchanges of particles, energy and in- 
formation with its surroundings. From an open quan- 
tum system point of view, the nonequilibrium electron 
dynamics is completely described by the master equa- 
tion through the physical observable {0{t)) = tr[Op(i)], 
where p[t) is the reduced density matrix of the cen- 
tral system that can fully depict the dynamics of elec- 
tron quantum states and is determined by the mas- 
ter equation. The transient electron transport should 
be in principle solved from the reduced density ma- 
trix as well. However, it has been struggled for many 
years without a very satisfactory answer in finding the 



2 



exact master equation for an arbitrary open quantum 
system. Most of the master equations used in the 
literature are obtained using semiclassical approxima- 
tion or perturbation truncation, such as the semiclas- 
sical Boltzmann equation^^'^^ or master equations un- 
der the Born-Markov approximation. ^^'^^ However, for 
the nanostrucutres with the extremely short length scales 
(~1 nm) and extremely fast time scales (~1 fs), the semi- 
classical Boltzmann equation or the master equation un- 
der the Born-Markov approximation are most unlikely 
applicable. 

Literarily, there are two fundamental but equivalent 
methods in dealing with modern nonequilibrium physics. 
They are the Schwinger-Keldysh nonequilibrium Green 
function tcchniquc'^''^''^ and the Fcynman- Vernon influ- 
ence functional approach. ^'^ In the practical applica- 
tions, both methods have their own advantages. The 
Schwinger-Keldysh nonequilibrium Green function tech- 
nique is a standard method for systematic perturba- 
tion calculations to various nonequilibrium physical pro- 
cesses in many-body systems. ^''■'^^ As it is evidenced 
the nonequilibrium Green function technique provides 
an c;xtremely powerful tool in the study of the steady- 
state transport properties in nanostructures.^'*''^^"^^ In 
the steady-state limit where the initial state of the nanos- 
tructurc is irrelevant, the transport current can be ex- 
pressed in terms of the nonequilibrium Green functions in 
the frequency domain which largely simplifies the trans- 
port problem. Without the time-dependence, the physics 
of the nonequilibrium transport is essentially determined 
by the density of states (the difference between the re- 
tarded and the advanced Green functions). The nonequi- 
librium Green function technique has been used to suc- 
cessfully describe a variety of new phenomena, such as 
Kondo effect, Fano resonance and Coulomb blockade ef- 
fects in quantum dots. The extension of the nonequilib- 
rium Green function techniques to the time-dependent 
transport phenomena in nanostructures has also been 
developed^°~^^'^^ where the time-dependent nonequilib- 
rium Green functions must be maintained. Except for the 
wide band limit (WBL),^^'^^ the relevant Green function 
calculations become rather complicated in order to tak- 
ing into account the different initial quantum state effect 
and the non-Markovian memory structure. ^'^'^'^ 

On the other hand, the Feynman- Vernon influence 

functional approach^'' are mainly used to study dissipa- 
tion dynamics in quantum tunneling problems''^ and de- 
coherence problems in quantum measurement theory. ^'^ 
It is very powerful to derive nonperturbatively the mas- 
ter equation for the reduced density matrix of an open 
system by integrating out completely the environmental 
degrees of freedom through the path integral approach, 
where the non-Markovian memory structure is mani- 
fested explicitly. In the early applications of the influence 
functional approach, the master equation was derived for 
some particular class of Ohmic (white-noise) environment 
for the quantum Brownian motion (QBM) modeled as 
a central harmonic oscillator linearly coupled to a set 



of harmonic oscillators simulating the thermal bath."*^' 
The exact master equation for the QBM with a general 
spectral density (color-noise environments) at arbitrary 
temperature that can fully address the non-Markovian 
dynamics was the Hu-Paz-Zhang master equation. Ap- 
plications of the QBM exact master equation cover var- 
ious topics, such as quantum decoherence, quantum-to- 
classical transition, quantum measurement theory, and 
quantum gravity and quantum cosmology, ctc.^'^'^'*'*^ 
Very recently, such an exact master equation is fur- 
ther extended to the systems of two entangled bosonic 
fields'*^"'** for the study of non-Markovian entanglement 
dynamics.*^ Nevertheless, using the influence functional 
approach to obtain the exact master equation has been 
largely focused on the bosonic (thermal) environments in 
the past half century. 

Meanwhile, using master equation to study quantum 
electron transport has also been attracting much atten- 
tions recently. With the nonequilibrium Green function 
technique, the master equation can be formally obtained 
in terms of real-time diagrammatic expansion up to all 
the orders, although practically most of the mas- 
ter equation approach used in quantum transport by 
far only take the perturbative theory up to the second 
order. On the other hand, with the help of the 
influence functional approach, an alternative formalism 
for the master equation for studying quantum transport 
and electron dissipative dynamics has recently been de- 
veloped in term of the hierarchical expansion. This new 
master equation approach provides a very efficient tool 
for the numerical study of the quantum transport prop- 
erties, including the accurate evaluations of the Coulomb 
blockage and Kondo transition dynamics. 

In this paper, we shall develop a nonequilibrium the- 
ory to describe nonperturbatively the transient electron 
transport dynamics intimately entangling with the elec- 
tron quantum decoherence in nanostructures. By extend- 
ing the influence functional approach to the fermionic 
reservoirs, we have recently derived an exact master 
equation for a large class of nanostructures.^^ This ex- 
act master equation is capable for the study of the full 
non-Markovian decoherence dynamics of electrons in na- 
noelectronics devices. In the present work, we shall ex- 
tend the previous work to the nanostructures in which the 
energy levels of the central region, the couplings to the 
leads and the external biases applied to leads are all time- 
dependent. We then derive the exact transient current 
in terms of the reduced density matrix within the same 
framework. The resulting transient current provides an 
exact non-linear response theory that can be used to in- 
vestigate the transient quantum transport processes ac- 
companied explicitly with the non-Markovian quantum 
relaxation and dephasing dynamics. The nonequilibrium 
steady-state transport theory based on the nonequilib- 
rium Green function technique that is widely used in the 
literature can be recovered as a long time limit from the 
present theory. 

The remainder of this paper is organized as follows. 
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In Sec. II, we outline the derivation of the exact master 
equation we obtained recently,^^ with the extension to 
the time-dependent Hamiltonian for a general nanostruc- 
ture. In Sec. Ill, we derive in detail the exact transient 
current in terms of the reduced density operator based 
on the Feynman- Vernon influence functional approach. 
The resulting transient current Ia{t) through the lead 
a is given by the general formula, Eq. (19), or equiva- 
lently Eq. (30) , which is valid for arbitrary bias and gate 
voltage pulses with arbitrary energy-dependent couplings 
between the central system and the leads. The relation 
between the reduced density matrix and the transient 
current is established through the master equation, see 
Eq. (23) where the superoperators £j(t) and £^ (t) given 
explicitly by Eq. (20) and (22) encompass all the back- 
reaction effects associated with the non-Markovian dy- 
namics of the central system interacting with the lead a. 
The steady-state current formula based on the nonequi- 
librium Green function technique can be recovered as a 
long time limit, see Sec. III.E. Applications to transient 
dynamics with Lorentzian-type spectral density and var- 
ious time-dependent ac bias voltages are demonstrated 
in Sec. IV. The summary and prospective are given in 
Sec. V. In Appendix, a close connection between the 
Feynman- Vernon influence functional approach and the 
Schwinger-Keldysh nonequilibrium Green function tech- 
nique for quantum transport theory is established. 

II. EXACT MASTER EQUATION FOR 
NANOSCALE ELECTRONIC DEVICES 

In this section, we shall outline the exact master equa- 
tion derived recently by two of us for nanoscale electronic 
devices^ ^ with the extension to the time-dependent en- 
ergy levels and couplings. We shall only list the necessary 
formula that we will use later for the derivation of the 
transient current. For more detailed derivation, please 
refer to Ref. [21]. We begin with the Hamiltonian of a 
prototypical mesoscopic system for electron transport, 

ij ak 

+ YiViM^)alcak + V*,,{t)cl^ai]. (1) 

iak 

Here the first summation is the electron Hamiltonian Hs 
for the central system where the electron-electron in- 
teraction is not considered. But as we have shown we 
can properly and explicitly exclude the doubly occupied 
states and the resulting master equation can be applied 
to the strong Coulomb repulsion (i.e. Coulomb block- 
age) regime. The second summation is the Hamiltonian 
J2a describing the noninteracting electron leads (the 
source and drain electrodes, ■ • ■ ) labeled by the index 
a (= L,R,---). The last term is the electron tunnel- 



ing Hamiltonian Ht between the leads and the central 
system. Since we will focus on the transient dynamics of 
quantum transport in this paper, the bias voltages Va ap- 
plied to leads are considered to be time-dependent. Thus 
the single electron energy levels in the leads is changed 
as Cak £ak{t) = Cak + (iVa{t)- Meanwhile, the en- 
ergy levels Cij of the central system and the couplings 
Viak between the central system and the leads are con- 
trollable through the gate voltages and external fields 
so that they can be, in general, also time-dependent: 
eij eij{t) = eij + Aij{t),Viak ^ ^iafc(t)- Through- 
out this work, we set ft = 1, except for the transient 
current where we will put ft back from its definition. 

The master equation for the nonequilibrium electron 
dynamics of the central system is derived based on the 
Feynman- Vernon influence functional approach. "'■^ As it 
is well known, the nonequilibrium electron dynamics 
of an open system is completely determined by the so- 
called reduced density matrix. The reduced density ma- 
trix is defined from the density matrix of the total sys- 
tem (the central system plus the leads, and the leads 
are treated as the reservoirs to the central system) by 
tracing over entirely the environmental degrees of free- 
dom: p{t) = trRptot(i)) where the total density matrix is 
formally given by ptot(i) = U{t,to)ptot{to)UHt,tQ) with 
the evolution operation U{t, to) = T cxp{—i j^^ H{T)dT}, 
and T is the time-ordering operator. As usual, we as- 
sume that the central system is uncorrelated with the 
reservoirs before the tunneling couplings are turned on:"^^ 
Ptot(io) = p(to) €5 PRito), and the reservoirs are initially 
in the equilibrium state: pR{to) = ie" /5«(^«-/i«Ar„) 
where f3a = l/{k^Ta) is the initial inverse temperature 
and Na = X)fc CafeCafe the particle number operator for 
the lead a. Then the reduced density matrix at an arbi- 
trary later time t can be expressed in terms of the coher- 
ent state representation^^'^" as 

xJ{^f,^'f,t\^o,^'o,to), (2) 

with I = (Ci,6r--) and i = (Ct,S=---) being the 
Grassmannian numbers and their complex conjugate de- 
fined through the fermion coherent states: ai\^) = 
and (^|a| = {^\^* ■ The propagating function in Eq. (2) is 
given in terms of Grassmannian number path integrals, 

J i'^']e^(^=[«' «-^*[«"'' ^'^^m-, i'a (3) 

where S'c[^,^] is the action of the central system in the 
fermion coherent state representation, and |'^'] is 

the influence functional obtained after integrating out all 
the environmental (reservoirs) degrees of freedom: 
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J to J kl 



+ E r'^^'(9ai,(r,rOC(rXK^O-9f«(T,T')(C(T)+CW)fe(r')+^K^')))}- (4) 

aij ■^*o 



The two-time correlations in Eq. (4) are defined by 
9aij{r,T') = $^Vi„fe(r)\A;,(r')e-/;"^-^-(-), (5a) 

k 
k 

(5b) 

which depict all the time-correlations of electrons in the 
leads through the interactions with the central system, 
and /a(eafe) = l/(e''«(*«'^~''«^ - 1) is the Fermi distribu- 
tion function of the lead a at the initial time to- 

This influence functional takes fully into account the 
back-reaction effects of the reservoirs to the central sys- 
tem. It modifies the original action of the system 
into an effective one, eA^'^^t^' I'C] = 
Q-kSenlii'i i ] which dramatically changes the dynamics 
of the central system. The detailed change is manifested 
through the generating function of Eq. (3) by carrying 
out the path_ integral with respect to the effective ac- 
tion Ses[$^; ^'^']. While the path integral V[$^; ^'^'] in- 
tegrates over all the forward paths ^{t),^(t) and the 
backward paths ^'(t),^'(t) in the Grassmannian space 
bounded by m = ^"/, ^(^o) = €o_and Cih) = i'o, ^'(t) = 
respectively. Since <S'eff[^4;^'^'] is only a quadratic 
function in terms of the integral variables, the path inte- 
grals in Eq. (3) can be reduced to Gaussian integrals so 
that we can use the stationary path method to exactly 
carry out these path integrals. The resulting generating 
function is:^^ 



^(4/>C/)^ICo,4o>*o) 



5^^|^exp{|,J,(i)4o 

+ ^fj2m'f + ^oJsmo + ^'oJlm'f}, (6) 



in which the time-dependent coefficients are given ex- 
plicitly as: Ji{t) = w'^{t)u{t), J2{t) = w^{t) - I, and 
J^{t) = u^{t)wl^{t)u{t) - J, with wf^{t) = [I- vl^{t)]-^. 
Here, we have also expressed the stationary paths in 
terms of N x N matrix variables u{t), vI^{t) and u{t) 
for to < T < t, where N is the total number of single 
particle energy levels in the central region, and the index 
/? implies that the corresponding function is tempera- 
ture dependent. They satisfy the following dissipation- 
fluctuation integrodifferential equations (i.e. the station- 



ary path equations of motion): 
u{t) + ie{T)u{T)+ ^ / dr'gair, t')u{t') = 0, (7a) 

U{t) + ie{T)u{T)- ^ / dT'g„(T, t')-u(t') = 0, (7b) 
v^{t) + ie{T)v^{T) +Y.r dr'g^iT, t')v^{t') 

a ''to 

= E f dT'gi{T,T')u{T'), (7c) 

a •'to 

with the boundary conditions u{to) — 1, u{t) = 1 and 
v^{to) = 0, where ga{T,T') and g^{T,T') are the non- 
local two-time correlation matrix functions, whose ma- 
trix elements are given by Eq. (5). As we will show 
in Appendix the matrix variables u{t), u{t) and v^{t) 
correspond respectively to the retarded, advanced and 
lesser Green functions, and ga{T,T') and ^^(t, r') are 
the retarded and lesser self-energies in the nonequilib- 
rium Green function technique. For the most general 
cases where all the parameters in the Hamiltonian (1) are 
time-dependent, the time translational invariance is bro- 
ken down. Then u{t) and tt(r) are usually independent 
except for the end-points where we have u{to) ~ M^(t). 
If all the parameters in the Hamiltonian of Eq. (1) are 
the time-independent, u(t) will become only a function 
of T — to and ■u(t) = u^t — t + to), as we have shown in 
Ref. [21]. 

Taking the time derivative of the reduced density ma- 
trix with the solution of the propagating function (6), 
together with the D-algebra of fermion creation and an- 
nihilation operators in the fermion coherent state repre- 
sentation, we arrive at the final form of the exact master 
equation we obtained previously^ ^ 



dpjt) 
dt 



= -i[H's{t),p{t)] 



+ 



(8) 



The first term (the commutator) in the master equation 

accoimts for the renormalized effect (including the time- 
dependent shifts of the energy levels and the changes of 
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the transition amplitudes between them) of the central 
system due to the interaction with the leads. The result- 
ing renormalized Hamiltonian is H'g{t) = J2ij 
While the rest terms in the master equation describe 
the dissipation and noise effects (which results in a 
non-unitary evolution of the central system) induced 
also by the interaction with the leads. All the time- 
dependent coefBcients in Eq. (8) are determined explic- 
itly by u{t),u^t) and v'^{t) as follows: 

a 

= ^E['*«(*)+'^«WJ'^-' (9b) 

a 

= Y.\-^i{t)+\i\t)],^. (9c) 

a 

Here we have used the relations obtained from Eq. (7), 
Ka{t)=f dTga{t,T)u{T)[uit)]-\ (10a) 

Jto 

Af(t) = fdT{g„it,T)vf'{r)-g0^{t,T)u{T)} 

Jta 

^ K,,{f)vf\t). (10b) 

As it is well known, the time-dependent coefficients 
in the master equation describe the non-Markovian dy- 
namics due to the interaction between the central sys- 
tem and the leads. The time-dependence of these coef- 
ficients is fully determined by the functions m(t),m(t) 
and v^{t) which are governed by the integrodifferential 
equations (7) where the time integrals involve the non- 
local time correlation functions of the reservoirs, Qaij, r') 
and g^{T,T'). These non-local time correlation func- 
tions characterize all the non-Markovian memory struc- 
tures of the central system interacting with its environ- 
ment through the coupling Hamiltonian Ht- By solving 
Eq. (7), one can completely describe the quantum deco- 
herence dynamics of electrons in the central region due 
to the entanglement between the central system and the 
leads. Thus, the master equation determines indeed the 
exact nonequilibrium dynamics of the system. It may 
be also worth mentioning that our master equation does 
not need to be in Lindblad form to preserve positivity 
since it is exact so that the dampening coefficients are 
time-dependent and the positivity is guaranteed. In the 
next section, we will derive the transient current passing 
through the central system within the same framework. 
The result will allow us to explore the transient dynamics 
of the nonequilibrium quantum transport accompanied 
with the non-Markovian decoherence phenomena, which 



receives more and more attentions recently due to rapid 
developments in nano-technology, spintronics and quan- 
tum information processing, etc.^'^ 

III. EXACT TRANSIENT CURRENT 

A. Influence functional approach 

The transient current from the Q-lead tunneling 
through a-j unction into the central region is defined in 
the Heisenberg picture as 

Ut) = -e ij^NUt)) = {{Mt),H{t)]). (11) 

Here, e is the electron charge, {0{t)) = Tr[0(t)p^t] '^ith 
tr = trgtrR, iV„(i) = Efc c^^fc Wcafc(i), and is the 
total density matrix in the Heisenberg picture. By ex- 
plicitly calculating the above commutation relation with 
the Hamiltonian of Eq. (1) and then transforming it into 
the Schordinger picture, we have 

Ut) = trs J2 - alAaiit)] , (12) 

i 

where the operators Aai{t) = tTn[J2k^ick{t)cakPtot{t)] 
and Al^-{t) = [Aai{t)]^ ■ The time-dependence of these 
two operators comes from the non-Markovian memory 
dynamics by tracing over the environmental degrees 
of freedom of the total density matrix ptot{t) in the 
Schrodinger picture. These two operators are indeed 
the effective (dressed) electronic creation and annihila- 
tion operators acting on the reduced density matrix of 
the central system, as a result of tracing over the reser- 
voir degrees of freedom for the corresponding operators 
acting on the lead a. 

Following the procedure of obtaining the reduced den- 
sity matrix through the influence functional approach,^^ 
we can write 

XJai{^f,^'f,t\^0,^0,to), (13) 

where the operator-associated propagating function is de- 
fined as 

(14) 

Similar to the calculation of the influence functional 

in Eq. (14), the operator-associated influence functional 
•^cdK^'f'^'] be calculated in the same way with the 
result, 

^ai m f = -i^ai [C, m^^; ^'^'] , (15) 
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where the functional expression of the effective electron 
annihilation operator is obtained as 

Aa[$,a = fdT{g^{t,r)${T)-g^„{t,T)[^{T)+^'{T)]}, 

(16) 

and the same influence functional is given by 

Eq. (4). It should be pointed out that the factorability 
of the operator-associated influence functional is due to 
the fact that the environmental path integral is exactly 
compiitable. 

Since the path integrals in Eq. (14) can also be reduced 
to Gaussian integrals, similar to the derivation of the 
master equation for the reduced density matrix (8), we 
use again the stationary phase method to exactly carry 
out the path integrals of Eq. (14). The result is: 

- iAai{$'f,^0, t)J{if,$'f,t\^o,$o, to), (17) 

where ^ai(^/,€o,i) = J^jiVavWioj + 2;mj(0^/j] with 
Vait) = JldTga{t,T)u{T) + Za{t)u{t), and Za{t) = 
[X'l{t) + Kait)v'^{t)][l-v'^{t)]-'^. The propagating func- 
tion J{^f,^'f,t\^o,^'Q,to) is given by Eq. (6) in the last 
section. Substituting Eq. (17) into Eq. (13), and using 
the identity, 

together with the D-algebra for the fermion creation and 
annihilation operators in the fermion coherent state rep- 
resentation, we obtain the eflfective electronic annihila- 
tion operator after tracing over completely the environ- 
mental degrees of freedom: 

Aai{t) = -lJ2 {-^aijWhpW + Pit)aj] + Kaij{t)ajp{t)}, 

(18) 

where the coeflicient matrices (t) and Ka{t) are the 
same as these appeared in the master equation (8) and 
are explicitly given by Eq. (10), p{t) is just the reduced 
density matrix determined by the master equation. 

Accordingly, substituting the above result into 
Eq. (12), the transient current can be directly calculated. 
The resulting transient current is rather simple: 

J„(i) = Tr{A^(t) + K^{t)p^'\t) + H.c.} 

= ~Rej\TTr[ga.{t,T)vf'{T)-g^{t,T)uiT) 

+ g^it,T)u{T)u-'{t)[p^'\t)-v^{t)]}, (19) 

where the notation Tr is the trace over the energy level 
basis of the central system. plj\t) = trs[a]aip(t))] is the 



single particle reduced density matrix. While the matrix 
elements oi ga{T,T') and g^{T,T') are defined by Eq. (5), 
and u{t), v^{t) and u{t) are determined nonperturba- 
tively by Eq. (7), as shown in the last section. Thus 
the transient current Ia{t) that flows from the lead a is 
completely determined within the framework of the mas- 
ter equation for the reduced density matrix. As we shall 
show later, from the above result we can easily reproduce 
the steady-state current at i ^ cxd in terms of the non- 
equilibrium Green functions and the Landauer-Biittiker 
formula. 



B. Relation between the transient current and the 
reduced density matrix 

In fact, the transient current defined by Eq. (12) can 
be written as a trace over the transient current opera- 
tor: Ia{t) = |trs[/c(i)], where the current operator is 
obtained directly from Eq. (18): 

Ut) = - E {^^^M^'UjPii) + alpit)a,] 

ij 

]^Ct{t)p{t). (20) 

Here we have introduced a current superoperator C'^{t) 
such that the current operator can be simply expressed as 
the current superoperator acting on the reduced density 
matrix. Note from Eq. (10) that Y.ai^» + ^a) = 7^ and 
Xla('*a + t^a) = 27, the transient current operator given 
by Eq. (20) is closely connected to the master equation 
(8) for the reduced density matrix p{t). 

Because of the trace over the states of the central sys- 
tem, the transient current of Eq. (12) can be alternatively 
expressed as 

Ia{t) = i^tT,Y^[aiAii{t) - Aa,i{t)al] = -|tr,l(t). 

(21) 

Then using Eq. (18) again, we have 

+ Kc,ij{t)ajp{t)4 + H.c.} = C-{t)p{t), (22) 

where C~{t) is defined as the superoperator for Ia{t). 
Now it is easy to check that the master equation (8) for 
the reduced density matrix and the transient current of 
(19) can be simply expressed as 

^ = -i[Hs{t),p{t)] + Y^[jrt{t) + C-itMt), 

(23a) 

I^t) = U^Ctmt)] = -U^C-{t)p{t)], (23b) 
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where Hg is the original Hamihonian of the central sys- 
tem. This analytical operator relation between the re- 
duced density matrix and the transient current shows 
explicitly the intimate connection between quantum de- 
coherence and quantum transport in nonequilibrium dy- 
namics. The reservoir-induced non-Markovian relaxation 
and dephasing in the transport processes are manifested 
through the supcropcrators of (20) and (22) acting on 
the reduced density matrix p{t). The time-dependent 
parameters X^{t) and Ka{t) appeared in the superoper- 
ators are given by Eq. (10) which are determined by the 
dissipation- fluctuation integrodifferential equation (7). 
This completes the nonequilibrium theory for transient 
electron dynamics in nanostructures we concerned. It 
should be pointed out that Eq. (23) was formally ex- 
pressed in terms of the real-time diagrammatic expan- 
sion and the hierarchical expansion. ^'^^ Here we have the 
analytical solution. 



The above equation reproduces exactly the transient cur- 
rent of Eq. (19). This also provides a self-consistent check 
for the expression of transient current derived directly 
from the influence functional approach. 

Taking the trace over the both sides of Eq. (25) and 
also noting the fact that N{t) = Tr p'^^\t), the total elec- 
tron occupation in the central system, we have 

^^-E^«W = -«), (28) 

a 

where /dis(i) is defined as the transient displacement 
current. It tells that sum over the currents flowing 
from all the leads into the central region equals to the 
change of the electron occupation in the central region, 
as we expected. In the steady-state limit t ^ oo, 
N{t) = Trp(i)(t) = so that /dis(t) = 0, as a conse- 
quence of current conservation. 



C. Current conservation law 

Now we should derive the current conservation law for 
a self-consistent check. Consider the equation of motion 
for the operator ata^ in the Heisenberg picture. 



\a-a 



^ejialai - eaalaj 



+ J2iVjc.k4c^k - CfcC^fca,). (24) 



ak 



Taking the expectation value of the above equation with 
respect to the state (the total density matrix in the 
Heisenberg picture) and then transforming it into the 
Schordinger picture, we obtain 

= m,p^'Ht)]+ij2ut). (25) 

a 

Here we have used the definition of the single-particle re- 
duced density matrix again, p^^\t) = tr[a]ai/9tot(i)] = 
tVs[a^jaip{t)], and also introduced a current matrix 

Ia]i{t) = itrs[Al^{t)aj - aJ^Qj(t)]. Equivalently, the 
transient current of Eq. (12) is simply given by Ia{t) = 
j-TTla{t). In other words, the equation of motion for 
p^^^{t) is directly related to the transient current. 

On the other hand, the equation of motion for the sin- 
gle particle reduced density matrix can also be obtained 
easily from the exact master equation (8). The result is 



dt 



uu 



-Vi)+p«(Mt)-i«t_y 



(1)1 - 



(1) 



H.c.) 



(26) 



Using the expression for 7^(t) given by (9c), and com- 
paring Eqs. (25) with (26), we have 



Ut) = - {A^(t) + Ko.it) p^'\t)+li.C.}. 



(27) 



D. Solution to single particle reduced density 
matrix and transient current 



Furthermore, we can rewrite Eq. (9c) as 



dt 



Comparing between Eq. (25) for the single particle re- 
duced density matrix p^^^{t) and the above equation for 
v^{t), it indicates that the solution of p^^\t) is just vl^{t), 
apart from an initial function. Explicitly, Eq. (25) and 
the above equation lead to 



dt 



It is easy to find from the above equation the following 
relationship between p^^\t) and v^{t): 

p^'Ht) = v^it) + u{t)p^'HtoWit), (29) 

where p^^\to) is the initial single particle reduced den- 
sity matrix. Accordingly, the transient current (19) is 
reduced to 

Ia{t) = ~Re J^dTTr[g4t,T)v^{T) - g^{t,T)uiT) 

+ g^{t,T)u{T)p^'\toW{t)y (30) 

The last term shows the explicit dependence on the ini- 
tial single particle reduced density matrix (including the 
initial electron occupation in each level and the electron 
quantum coherence between different levels in the cen- 
tral region). This term is an important ingredient in the 
study of transient dynamics for practical manipulation of 
a quantum device in real time and it will usually vanish 
in the steady-state limit t ^ oo. 
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E. Steady-state limit 

To make an explicit comparison with the steady-state 

current in terms of the non-equilibrium Green functions 
and the Landauer-Biittiker formula used in the litera- 
ture, we introduce the spectral density of the lead a: 
ray(t^) = '2TTj2k'^aikV*j,^5{uj - e„fe) and take a time- 
independent bias voltage explicitly. Then the two-time 
correlation functions of the a-lead can be written as 



, f cLu 

9a{T-T) = / —Ta{oj)e- 



■iu!(r—r') 



2^ 



-i(jj{r—r') 



(31a) 
(31b) 



where fa{i^) = l/(e''°^'^~^°'' + l) is the Fermi distribution 
function of the a-lead at the initial time to. Using the 
Laplace transformation, i.e., f{z) = dte~^^*~*^°^ f{t) 
with z = —iui, we have 



9a (<^) 



r duj' 

J 2^ 



r«(a;') 



iS^(a;), (32) 



i.e., S'a(a;) is the retarded self-energy induced from the 
lead a. The Laplace transformation of Eq. (7a) for u{t) 
gives 



S'-(a;) 



(33) 



where G^{w) is just the retarded Green function, and 
sums over S'j(ti') for all a. The advanced Green 
function is simply given by u{oj) = — iG"(a;) = m^(cl)). 
Furthermore, for the time-independent Hamiltonian, the 
explicit solution of Eq. (7c) is 

V^{t)= f dn f dT2U{Ti)g'^{T2-n)u\T2). (34) 

Taking the long time limit: t oo, its Laplace transfor- 
mation gives 



= -iG''(a;)S<(a;)G"(a;) 



-iG<{uj). (35) 



Here we have also used the relation g^{(-o) = — iS<(ti>). 
More explicit relations between u{t).u{t) and v^^{t) with 
the retarded, advanced and lesser Green functions in the 
real-time domain are presented in Appendix. 

Substituting above results into Eq. (30) in the steady- 
state limit t ^ cx), we obtain the steady-state single par- 
ticle reduced density matrix and current: 

P^^= J ^G'^(^^)[X^r„(a;)/„(a;)]G»(a;), (36a) 

7„,st = ylm J gTv{r„(a;)[G<(a;)-/«(u;)G»(a;)]}. 

(36b) 



This reproduces the steady-state current in terms of the 
nonequilibrium Green functions in the frequency domain 
that has been widely used. If we consider specifically a 
system coupled with left (source) and right (drain) elec- 
trodes, i.e. a = L and R, respectively, and also assume 
that the spectral densities for the left and right leads have 
the same energy dependence: rz,(a;) = XTji{u)), where A 
is a constant. Then the net steady-state current flowing 
from the left to the right lead is given by 

^.. = 1/ ^[/tM-/;,M]r(„). 

This is the generalized Landauer-Biittiker formula.^''' 
Thus, we have simply recovered the nonequilibrium 
transport theory at the steady-state limit. In fact, we can 
also reproduce the transient transport theory in terms of 
the nonequilibrium Green function technique, as shown 
in Appendix. 

We now summarize the main results derived in this 
Section. We obtain the general formula of the tran- 
sient current Ia{t) through the lead a which is given 
by Eq. (19) or equivalently Eq. (30), accompanied with 
the exact master equation for the reduced density ma- 
trix p{t). This general transient current is valid for ar- 
bitrary bias and gate voltage pulses with arbitrary cou- 
plings between the central system and the leads. We 
also establish explicitly the connection of the transient 
current with the reduced density matrix, i.e. Eq. (23), 
through the superoperators £+ (t) and C~ (t) determined 
by Eq. (20) and (22), which encompass all the back- 
reaction efFe(;ts associated with the non-Markovian dy- 
namics of the central system interacting with the lead a. 
The current conservation law is also explicitly derived. 
The steady-state current formula based on the nonequi- 
librium Green function technique is reproduced as a long 
time limit. These results allow us to explicitly analyze 
the time-dependent quantum transport phenomena inti- 
mately entangling with the electron quantum coherence 
and non-Markovian dynamics through the reduced den- 
sity matrix. The latter describes completely the evolu- 
tion of electron quantum coherence inside the nanoelec- 
tronics device. Therefore, once one solves the equation 
of motion (7) for u{t),u{t) and v^{t), the full nonequi- 
librium dynamics of the nanostructures can be analyzed 
explicitly. 



IV. ANALYTICAL AND NUMERICAL 
ILLUSTRATIONS 

A. Time-independent bias voltage in the WBL: an 
analytic solution for transient dynamics 

For practical applications, we take a simple system as 
an illustration: a single dot containing only one spinless 
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level that couples to the left and right leads. The Hamil- 
tonian of the dot is simply written as H = ea^a. This sys- 
tem only contains two states, the empty state |0} and the 
occupied state |1). All the corresponding matrixes (de- 
noted by the bold symbols) in the above formula are then 
reduced to a single function, such as u{t) = u{t), v^{t) = 
vl^it), andp(i)(t) = = pii(i) = N{t). We will first 

consider a time-independent bias voltage V. For simplic- 
ity, we also assume that the tunneling couplings between 
the leads and the dot as well as the densities of states for 
the leads are energy-independent. In other words, the 
spectral density becomes a constant Ta{Lo) = Ta- The 
non-local time correlation functions is reduced to 



9a{T - t') = Va5{T - t'), 



(38a) 
(38b) 



This corresponds to the wide band limit (WBL) in the 
literature. 

To be explicit, we take the initial time to = and let 
e = ft = 1 in the following calculations. In the WBL, the 
solution of Eq. (7) is 

u{t) = exp I - (ie + ^)r| , u{t) = u^t - r) (39a) 

27r (e - c^)2 + (r/2)2 
X {e-"-2e-"/2cos[(e-a;)t]}, (39b) 



where F = -|- F/j and is the solution of v^{t) at the 
steady-state limit: 



^LfLjm) + TRfii{u)) _ (1) 

27r (e-t^)2 + (F/2)2 ~ ' 



(40) 



The electron occupation of the dot is calculated by 
Eq. (29) as 



N{t) = p^^\t) = e-^*p^'\(])+v^\t), 



(41) 



where A^(0) = ^^^•'(0) is the initial electron occupation 
of in the dot. Obviously, in the steady hmit, Nst = w^. 
This is not surprised since p^^^ (t) and (t) obey the same 
equation of motion that must lead to the same result in 
the steady-state limit. 

The transient current can then be analytically ob- 
tained: 



Iait)=Ia.,st-Ta[N{t)-N,t] 



- e 2 



r|Fcos[(e-w)i] 



27r (e - + (F/2)2 

-2(e-a;)sin[(e-a;)t]|, (42) 

where the steady-state current is 
J 2n 



(e-a;)2 + (F/2)2 



Due to charge conservation, it is necessary to check the 
displacement current which obeys the relation (28). In 
the WBL, 



Jdis(t) = Fe 



doj TLfL{uj)+TRfR{u;) 



(e-a;)2 + (F/2)2 
X { - Fe-" + Fe-"/2 cos[(e - u)t] 

+ 2(e-a;)e-"/2sin[(e-a;)t]} 
dN{t) 



dt 



(44) 



At the steady-state limit, the displacement current 
Idis{t ^ oo) = 0, as a consequence of current conser- 
vation. It is also straightforward to calculate the net 
current: 

Inetit) = hit) - lR{t) 

= I,,-{TL-Tn)[Nit)-Nst] 
_ f doj TL.fL{i^) - TRfji{u}) 



I 



27r (e - a;)2 + (r/2)2 

X |f cos[(e - uj)t] - 2(e - u) sin[(e - (j)t] | 

(45) 



where the stationary net current is 



I,t = 2TlTr 



I 



277 (e - w)2 + (F/2)2 



(46a) 



As we can see, once we solve the dissipation-fluctuation 
integrodifferential equations, Eq. (7), a complete time- 
dependence of all the physical quantities, such as the 
electron occupations in the central system and the cur- 
rents flowing from each lead to the central region, can be 
obtained with the explicit dependence of the time and 
the initial electron occupation without ambiguity. Prom 
Eq. (42), we sec that the initial current /„(0) = -r„iV(0) 
which depends on the initial occupation of the dot. This 
result is consistent with the electron occupation in the 
dot, Eq. (41). For zero initial occupation, the initial cur- 
rent is zero. Note that some of the above results are also 
obtained recently using the nonequilibrium Green func- 
tion technique2°. 



B. Non-Markovian memory structure 

Realistically, the spectral density of the leads must de- 
pend on the energy. Here we take the energy dependence 
as a Lorentzian-type form:^^'^^'2i 



F W'^ 



(47) 



(43) 



where F„ describes the coupling strength and is the 

line width of the source (drain) reservoir with a = L{R). 
Obviously the WBL, Fc(a;) = F^,, is achieved by simply 



10 



letting Wa oo. The lead correlation functions with 
time-independent voltage can be parameterized as^^ 



M 



(48a) 
(48b) 



The first term in Eq. (48b) with m ~ arises from the 
pole of the spectral density function, with 



1 + e-^0c,w^ 



7aO = Wa+ i^la■ (49) 



The other terms with m > (M — > oo in principle) arise 
from the Matsubara poles, where the relevant parameters 
are explicitly given as 



riam = —Ta{-i^am), m = l,---oo, (50a) 

Ha 



_ (2m- l)7r 

Pa 



(50b) 



There are four typical timescales in an open system to 
characterize the non-Markovian memory structure: the 
timescale of the central system ('^ 1/e), the timescale 
of the reservoirs 1/VFq), the mutual timescale due to 
the coupling between the system and the reservoir ['^ 
I/Fq), and the thermal timescale (~ (3a or 1/^q). The 
timescale of the central system is usually fixed when the 
system is set up. Then comparing with the character 
time of the central system, a smaller finite line width Wa , 
a relatively large coupling Fq,, a low temperature 1//3q, 
and a comparable bias voltage eV = fJ-L — fJ'R to the 
energy levels of the central system will lead to a stronger 
non-Markovian memory processes for a Lorentzian-type 
spectral density, as we have demonstrated in [21]. 

In fact, the line width Wa in a Lorentizan-type spectral 
density is the main factor leading to the non-Markovian 
dynamics in transient transport. In the WBL, Wa —>■ oo, 
the dominated memory structure are mostly washed out. 
This can be seen directly from the reservoir correlation 
functions. The correlation function ga{t — T) of Eq. (48a) 
can be simplified to ^6{t — t) in the WBL. While for 
g^{t — t) in Eq. (48b), the first term (m = 0) is also 
simplified to a delta function oit — r but the other terms 
(to > 1) are apparently changed not too much: 



9i{t-r) 



Fa 



2(l + e-*/3=>^=) 



5{t-r) 



M 



,-7=,r.(t-r) 



(3a 



(51) 



The profile of this temperature-dependent time correla- 
tion function is plotted in Fig. 1. If we take further a 
high temperature limit /Sq — > 0, the summation term 
in Eq. (51) will also be reduced to a delta function of 
t — r, see Fig. 1. Then no any memory effect remains. 
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FIG. 1: (Color online) 3D plot of temperature-dependent non- 
local time correlation function g1{t — r) (in unit: F^) as a 
function of temperature for hl ~ 2.5T and Tl ~ 0.5T. 
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FIG. 2: The temperature-dependent non-local time correla- 
tion function {t — r) with different bias voltages at a fixed 
temperature /3"^ = k^T = O.IF and Fl = 0.5F. The plots 
of (a) and (b) for fj,L ~ 0.5F, (c) and (d) for fiL = 30F (the 
large bias limit). 



and a true Markov limit is reached at high temperature 
limit. On the other hand, for a large bias voltage limit 
eV = /iL ~ I^R — > oo, we have /L(t^ + 3-) ~> 1 and 
fni^ — ^) which lead to g^{t — r) ^ ^^(t — r) = 
^S{t - r) and g'^{t - r) ^ in the WBL. This re- 
duces the electron dynamics into a Markov limit again, 
as it has been widely used in the literature. ^^"^^ How- 
ever, for a relatively low temperature or a finite bias 
voltage, there appears to exist some non-Markovian ef- 
fect in the WBL, coming from the summation term in 
Eq. (51). Fig. 2 shows the time dependence of the cor- 
relation dependence with different voltages. As one can 
see, the temperature-dependent time-correlation function 
does not approach to a delta function in time. 

To examine if some non-Markovian effect can still sur- 
vive in the WBL, we numerically calculate the transient 
current passing through the single level resonant tunnel- 
ing nanostructure considered in the last subsection. In 
the sequential tunneling regime /xl > e > ijlr, the exact 
solutions to the occupation and the transient current are 
close to the Markov limit (differing by a few present ex- 
cept for a very short time scale from the beginning), as 
shown in Fig. 3(a)- (d). While, in the co-tunneling regime 
with fXa S> e, the exact solution of the occupation and 
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FIG. 3: (Color online) Comparison of the non-Markovian dy- 
namics with Markov limit in the WBL for the electron occu- 
pation (left row) and the transient current (right row). Here 
we take the parameters Tl = — O.SF, f3~^ = ksT — O.IF, 
HL.R = ±eV/2 = ±5F, and (a)-(b): eV = lOF, e = 2F; 
(c)-(d): eV = 20F, e = 2F; (e)-(f): eV = lOF, e = -lOF. 



the current are still almost the same as their Markov 
solution except for the very short time scale from the be- 
ginning, see Fig. 3(e)-(f). These results indicate that the 
WBL (an extremely short character time of the reser- 
voirs) will dominatively suppress the thermal timescales. 
In other words, when the line width W —>■ oo, not only 
for the high temperature and large bias limit, but even 
for a finite temperature and a finite bias voltage, the 
non-Markovian effects becomes quite weak that are most 
likely negligible in experiments. Thus the WBL mainly 
takes into account the Markov dynamics. The manifesta- 
tion of the non-Markovian memory structure then should 
go beyond the WBL.^^ In Fig. 4, we calculate the tran- 
sient current with different line widths to demonstrate 
the non-Markovian effect in transport phenomena. As 
we can see, the transport dynamics is significantly differ- 
ent from the Markov limit for a small W. Increasing the 
value of W will decrease the memory effect accordingly. 
When the line width W > 50r, the exact solution of the 
transient current closely approaches to the Markov result 
that is consistent with the WBL. A analysis of the non- 
Markovian dynamics in this simple resonant tunneling 
system has also recently been studied using Heisenberg 
equations of motion. 




2 4 6 
Time t (l/F) 

FIG. 4: (Color online) The transient current Jh^Rit) through 
the left and right leads with a Lorentzian-type spectral density 
for different bandwidths: dashed line, the WBL (W = oo), 
other lines, from the top to bottom for (a) and from the bot- 
tom to top for (b), correspond to W = 50F, 20F, lOF, 5F, 2F, 
respectively. Other parameters we used are Fl = F_r = O.SF, 
= k^T = O.IF, ^iL,R = ±eV/2 = ±5F, and e = 2F. 



transport dynamics is also a central ingredient in many 
different experiments, such as single-electron pumps and 
turnstiles with time-modified gate signals moving elec- 
trons one by one through quantum dots,^''^^^ and the 
study of the quantum capacitance and inductances with 
ac voltage response. ^^"^^ All these problems can be stud- 
ied explicitly in the present theory now. The correspond- 
ing nanoscale devices can be modelled as a resonant tun- 
neling nanostructure considered in this section. The ap- 
plying time-dependent voltages are taken to be the most 
commonly interesting ones. For time-dependent voltages, 
as we mentioned in Sec. II the single particle energy lev- 
els of the leads are changed to eak{t) — £ak + eVa{t). The 
non-local time correlation functions of the leads can be 
expressed as 



ga{T,T') ^ exp^^-ie J driVa(ri)| ^^(t - r'), (52a) 
5f(r,r') = exp|-ze^^ driy„(ri)| (r - r'). (52b) 



C. Transient transport dynamics with 
time-dependent bias voltage 

We are now in the position to study the transient trans- 
port phenomena in response to time-dependent bias volt- 
ages. The calculation of the time-dependent electron cur- 
rent in response to time-dependent bias voltages has re- 
ceived a great attention recently, through various differ- 
ent theoretical approaches. Transient 



With the parametrization of Eq. (48), it is not difhcult to 
numerically calculate the transient electron dynamics for 
arbitrary line width Wa- In the following calculation, we 
set the symmetric ac voltages, i.e., //^(t) = eV{t)/2 and 
/ifl(t) — —eV{t)/2. For the quantum dot with a single 
level, the reduced density matrix can be fully character- 
ized by the electron occupation in the dot. The exact nu- 
merical results for the transient current and occupation 
due to different types of applied ac voltages are presented 
as follows. 
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FIG. 5: (Color online) The occupation and transient cur- 
rent response to exponential time-dependent voltage with r — 
1.5/r. The other parameters are eV = lOF, Tz, = = l.ST, 
f3-^ = ksT = O.ir, Wl = Wr = W^ 20r and e = 2r. 



FIG. 6: (Color online) The occupation and transient current 
response to oscillating voltage with the dc component eVo = 
lOr, the ac component eVc = 2r and the oscillating frequency 
ujc ~ 4r. The other parameters are the same as in Fig. 5. 



Exponentially time-dependent bias voltage 

We shall first study the transient transport dynamics 
in response to an exponentially time-dependent bias volt- 
age V{t) = V{1 — e"*/"^), where t > is a time dominat- 
ing switch-on rate of the voltage. The asymptotic limit 
r — > 0"*" corresponds to a step function. The numerical 
results are plotted in Fig. 5. The first peak shown in the 
displacement current arises from the cotunneling process 
during the initially short time scale. At beginning, the 
Fermi surface of the both leads are nearly equivalent to 
zero and the dot energy level is higher than the Fermi sur- 
face. The initial currents through both leads are equal to 
each other due to the totally symmetric structure. This 
leads to a zero initial net current. This feature is com- 
mon for other type of ac bias voltages discussing later. In 
general, the emergence of the peaks in the current corre- 
sponds to a steep transient behavior of the electron states 
(i.e. the occupation for the single-level dot) occurs inside 
the dot. The non- linear response to the time-dependent 
bias is clearly manifested in the change of both the elec- 
tron state in the dot and the transient current through 
the leads (including the individual current through each 
lead and the displacement and net currents), as we plot- 
ted in the Fig. 5. In particular, the net current changes 
in time that closely follows the change of electron occu- 
pation in the dot, while the displacement current depicts 
the steep change rate of the occupation, as we expected 
from Eq. (28). 

Oscillating bias voltage 

Now let us move to the transient dynamics driven by 
an oscillating voltage:^'* V{t) = Vq — VcCOs{uJct), where 
1^ is a dc component, and Vc and ujc are the oscillation 
amplitude and frequency of the ac component. The cor- 
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FIG. 7: (Color online) The transient current (black line) re- 
sponse to the oscillating voltage (the dashed red line) with 
eVo ~ eVc = lOr and the oscillating frequency uJc ~ IF for 
(a) and uic = O.SF for (b). The other parameters are the same 
as in Fig. 6. 



responding exact numerical solution is shown in Fig. 6. 
As one can see all the quantities, the electron state in 
the dot and the transient currents, have the similar os- 
cillation behavior as the ac voltage oscillation, where the 
steady-state values and the oscillation around the steady- 
state values are determined by the dc voltage Vq. The 
oscillation amplitude around the steady-state is propor- 
tional to the ac voltage amplitude Vc and the oscillation 
periodic is mainly given by T = l-KjiOc- However, the 
transient dynamics of the occupation and current does 
not always strictly follow the ac voltage oscillation. The 
current oscillation is a little slantwise compared to the ac 
voltage oscillation. With increasing the amplitude of the 
ac voltage Vc and decreasing the oscillation frequency Wc, 
a sideband oscillation occurs in the electron occupation 
as well as in the transient currents as shown in Fig. 7. 
Physically this sideband oscillation is induced by the si- 
nusoidal behavior of the ac voltage. It can be understood 
analytically under the WBL:^'' 
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FIG. 8: (Color online) The occupation and transient current 
responses to Gaussian voltage pulse with n = 3/r and T2 = 
0.5/r. The other parameters are the same as in Fig. 5. 



where the Bessel function satisfies J_„(z) = (— l)"J„(z) 
and Al^r = T^- 

Gaussian pulse 

The last example we shall study is the tran- 
sient dynamics droven by a Gaussian pulse V(t) = 

l^exp{ — ■'^ — T^}- The width and the center of the pulse 
are determined by T2 and ri, respectively. The exact 
numerical result plotted in Fig. 8 shows that the net cur- 
rent peak emerges (a slightly delay) just after the voltage 
pulse while the corresponding response of the electron 
occupation delays significantly. This behavior is easy 
to be understood that the external voltage pulse leads 
to a sharp change of the electron occupation due to the 
delay response. The shape change of the transient cur- 
rent comes from the largest change rate of the electron 
occupation in the dot. The delay response effect is de- 
termined by the tunneling rate F. Outside the voltage 
pulse, the tunneling current comes completely from the 
cotunneling effect and the occupation in the dot decays 
to a stationary value. It is interesting to note that the 
response of the occupation and the currents to the co- 
tunneling process is different before and after the voltage 



pulse. Before the voltage pulse, the system is dominated 
by the cotunneling process because the Fermi surface of 
the both leads are nearly equivalent but is lower than 
the dot energy level (hl — fJ-R — < ^ — 2F). The 
aligned fermi surfaces double the peak amplitude of the 
displacement current which also gives the zero net cur- 
rent, while the occupation has a corresponding change 
regarding the change of the displacement current. With 
the voltage increasing, the Fermi surface of the left lead 
moves up over the dot energy level while that of the right 
lead moves down below the dot energy level such that 
the system gradually approaches the sequential tunnel- 
ing regime. This drives the electron flowing from the left 
lead to the dot and then to the right lead. Such a process 
results in a sensitive response of all physical quantities, 
the electron occupation in the dot, the displacement and 
the net currents. Then with the voltage decaying to zero, 
the electron residing in the dot is favorable to cotunnel 
to the left lead which gives a negative current and finally 
reaches the steady-state. The transient electron dynam- 
ics with the non-linear response to this Gaussian pulse is 
in particular useful for quantum feedback control to the 
electron states in the dot through the transient current 
that we will study in the future. It is also favorable to 
study the quantum capacitance and inductances. ^^'^^ 

Note that the above transient dynamics starts with a 
zero initial occupation, i.e., /Ooo(^o) = 1 and pii(to) = 
1 - Poo(io) = [N{to) = 0]. This implies that the last 
term of Eq. (30) which is often ignored in the noncqui- 
librium Green function technique^ ^'^^ has no contribu- 
tion to the transient current in the above numerical cal- 
culations. If the dot is initially occupied, the transient 
current will be quite different although the steady-state 
limit is the same because the initial electron distribution 
vanishes in the steady-state limit. In Fig. 9, we show 
the exact numerical result with such a situation where 
we take the simple step-pulse voltage as an example. It 
shows that the initial occupation in the central region 
has measurable contributions in studying the transient 
dynamics, especially for the ultrafast (extremely short 
time) operations in a quantum device. 

Combining all these analysis together, as one can see 
the exact numerical solutions presented here have demon- 
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FIG. 9: (Color online) The occupation and transient current 
response to step pulse voltage (i.e. a constant bias after t = 0, 
the dot-dashed line) for the initial condition of pii(to) ~ 
(solid lines) and pii(to) ~ 1-0 (dashed lines), respectively. 
The other parameters are the same as in Fig. 5. 



strated that both the electron states in the dot and the 
transient currents past through it have a clear non-hnear 
response to the external bias voltage pulses, in partic- 
ular within a short time scale after the pulse is turned 
on. These ultrafast non-linear response properties will 
provide very useful information for the manipulation of 
the device states as well as the quantum feedback con- 
trols for practical applications. However, since we only 
consider here a very simple device with a single-level dot, 
the quantum coherence and decoherence dynamics to the 
transient current do not be manifested in these numerical 
solutions. To demonstrate the quantum coherence prop- 
erties in the transient transport dynamics, it is necessary 
to have the device containing at least two levels (or a 
single level with electron spin degrees of freedom) in the 
central region. Also, in the above numerical calculations, 
we take a rather large finite line width, Wa = 20r. The 
non-Markovian memory structure will become more sig- 
nificant when the line width becomes narrower. ^^'^^ We 
will examine in detail these features within the present 
theory in our future work. 



V. SUMMARY AND PROSPECTIVE 

In summary, we have established a nonequilibrium the- 
ory for the transient quantum transport dynamics via 
the Feynman- Vernon influence functional approach. We 
obtain an analytical relation between the master equa- 
tion, Eq. (8), and the transient current, Eq. (19), de- 
termined by the same time-dependent (non-Markovian) 
coefficients. In particular, the master equation and the 
transient current are explicitly related to each other in 
terms of the superoperators acting on the reduced den- 
sity matrix of the central system, see Eq. (23). The 
back-reaction effect of the gating electrodes to the cen- 
tral system is fully taken into account by the time- 



dependent coefficients in the master equation (8) through 
the dissipation-fluctuation integrodifferential equation 
(7). The non-Markovian memory structure is non- 
perturbatively built into the integral kernels in these 
equations of motion. The resulting transient current, 
Eq. (19) or (30), is rather simple and it recovers the 
steady-state current in the nonequilibirum Green func- 
tion technique and the Landauer-Biittiker formula at the 
long time limit. This exact nonequilibrium formalism 
should provide a very intuitive picture of how the change 
of the electron quantum coherence in the central system 
is intimately related with the electron tunneling processes 
through the leads, and therefore non- linearly responses 
to the corresponding external bias controls. This the- 
ory is applicable to study a variety of quantum transport 
phenomena involving explicitly non-Markovian quantum 
decoherence behaviors, in both stationary and transient 
scenarios, at arbitrary temperatures of the different con- 
tacts in the weak Coulomb interaction regime. 

As a simple illustration, we apply the theory to a 
simple model of electron tunneling though a single level 
quantum dot. We obtain all analytical solutions in the 
wide band limit (WBL) where the dependence of the 
initial electron occupation in the dot is explicitly man- 
ifested. Taking a more realistic spectral density with 
Lorentzian shape, we show that the Markov limit is a 
good approximation in the WBL. The non-Markovian 
memory effect is dominated by a finite line width of the 
spectral density. Under the ac bias voltage pulses (in- 
cluding the step pulse, the Gaussian pulse and the os- 
cillation pulse), we have demonstrated the ultrafast non- 
linear response of the electron occupation and currents 
to the ac bias. We find that the current evolutions are 
more sensitive to the energetic configuration of the dot 
than the occupation evolution. This feature is very useful 
for quantum feedback controls in quantum information 
processing. More applications will be presented in the 
future works. These include the decoherence transport 
dynamics in quantum dot devices such as quantum dot 
Aharonov-Bohm inteferometers; the nonequilibrium dy- 
namics and the real time monitoring of spin polarization 
processes in nanostrucutres; also the transient transport 
dynamics in molecular electronics, as well as the appli- 
cation to bio-electronics such as DNA junctions, etc. 

In the end, we should also point out that the present 
theory is developed without considering the electron- 
electron interaction in the central region and therefore it 
is mainly valid in the weak Coulomb interaction regime. 
It is not difficult to extend the present theory to the 
strong Coulomb Blockage regime by properly excluding 
the doubly occupied states in the central region of the 
nanostructure, as we have shown explicitly in [21]. Al- 
though studying the above two extreme limits, the ex- 
tremely weak and the extremely strong Coulomb inter- 
action regimes together, could reveal a significant under- 
standing to various quantum transport phenomena, there 
is increasing discussion for the intermediate Coulomb in- 
teraction regime^^'^^ where the analysis of transport dy- 



15 



naniics should become much more comphcatcd. In this 
situation, the path integrals of Eqs. (3) and (14) may 
not be carries out exactly, and therefore it is not easy 
to find an exact master equation and an exact expres- 
sion of the transient current. But the master equation, 
Eq. (8), and the transient current, Eq. (19), can still 
serve as a good approximation with respect to the sad- 
dle point approximation or loop expansion, where the 
Coulomb interaction must be included sclf-consistcntly 
in the dissipation-fluctuation integrodifferential equation 
(7). This approximate treatment could provide a basic 
understanding to the quantum transport phenomena in 
the intermediate Coulomb interaction regime. The de- 
tailed extension of the present theory to the interacting 
electron systems is in progress now. Nevertheless, the an- 
alytical nonequilibrium theory we presented in this paper 
has the advantage of combining the decoherence dynam- 
ics with transient transport dynamics together to explore 
time-dependent physical phenomena and may also pro- 
vide a basic theory for quantum feedback controls in var- 
ious nanoelectronics devices. 
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APPENDIX: THE RELATIONS BETWEEN 
u{t),u{T), v'^it) AND THE NONEQUILIBRIUM 
GREEN FUNCTIONS 

As we see both the master equation (8) and the tran- 
sient current (30) are completely determined by the prop- 
agating matrices of the stationary paths: m(t), u{t) and 
v^{t). Here we shall show that these propagating ma- 
trices are directly related to the retarded, advanced and 
lesser Green functions in the nonequilibrium Green func- 
tion technique. In our previous work,^^ the propagating 
matrices u{t),u{t) and v^{t) were introduced to sim- 
plify the stationary path equations of motion (with the 
convention x = yz for Xi = J2j Uij^j) ^ follows: 

^(r) = «(r)$(to) + v^irMt) + ^'(0], (A.la) 
m+^'{r)=u{TMt)+^'{t)]. (A.lb) 

In fact, Eq. (A.l) shows that u{t) is a propagating ma- 
trix of the forward stationary paths ^(t) starting at to, 
while v^{t) mixes the forward path ^(t) and the back- 
ward path $'{t) started backwardly from t, and u{t) is 
a backward propagating matrix of the stationary paths. 
In fact, Eq. (A.l) shows that the transformation matrices 



u{t), u{t) and v^{t) should be defined more precisely as 
u{t) = u(T,to), u{t) = ut(t,r) and t;^^(r) = v'^^(r,f): 

^(t) = u(r,to)^(to) + v'^(r,t)K(f) +^'(t)], (A.2a) 
|(r)+^'(r) = ut(t,r)K(t)+^'(t)], (A.2b) 

where the Grassmannian variables ^(r) and ^'{t) rep- 
resent the forward and backward electron paths in the 
functional path integrals. Then Eq. (A.l) directly tells 
that u(t, to) describes the electron propagation (repre- 
sented by ^(r) in the Grassmannian space) from the ini- 
tial time to to the time r so that it is just the retarded 
Green function, namely, 

u{T) = u{T,to)=lG^{T,to) 

= 0{T-to){{ai{T),al{to)}). (A.3) 

Eq. (33) is a justification for this relation. While, u^(f,r) 
describes the inverse propagation of the electron [or the 
backward propagation represented by ^'(r)] from the 
time t to the time r such that it is indeed the advanced 
Green function: 

u{t) ^ u\t, t) = -iCir, t). (A.4) 

Meantime, Eq. (7a) indicates that the time correlation 
function of the a-reservoir: 

gaij{Tl,T2) = iS^y(ri,r2) 
k 

(A.5) 

as a back-reaction effect of the a-lcad to the central sys- 
tem, is the retarded self-energy. These relations are also 
justified by Eqs. (7a- 7b). 

The function \/^{T,t) describes the electron propaga- 
tion mixing the forward and backward paths so that 
it is related to the lesser Green function defined by 
^i)(^'^) ^ *('^j(i)0'i(''')) in the nonequilibrium Green 
function formalism. In fact, it is not difficult to find 
the explicit solution of Eq. (7c): 

t,/3(r)=V^(r,t)= Tdn f dT2U{T,n)g^{n,T2Wit,T2) 

J to J to 

(A.6) 

which in terms of Green functions becomes 

sj^{T,t)=-i j dn f dr2 G''(r,ri)S<(Ti,r2)G'^(T2,t), 

'^to J to 

(A.7) 

where 

S<(ri,r2) = i5^(ri,r2), (A.8) 

is the lesser component of the self-energy. On the other 
hand, the single-particle reduced density matrix is related 
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to the lesser Green function by p'^^\t) = —iG^{T^t)\r=t- 
From the relation of Eq. (29), we find that v^{t) can be 
expressed in terms of the lesser Green function as follows: 

G<{T,t) =l[u{T)p^^\toW{t)+V^{T)] 

+ r dn I dT2 G'-(T,Ti)S<(ri,r2)G'^(r2,i). 

J to J to 

(A.9) 

This provides indeed a general solution of the lesser 
Green function in the nonequilibrium Green function 

technique. In the previous investigation of transient dy- 
namics, the first term in Eq. (A.9) that sensitively de- 
pends on the initial electron distribution in the central 
regions are often ignored. '^-'^''^^'^^ In fact, the second term 
in Eq. (A.9) can only be identified as the lesser Green 
function G<(t, t) in the steady-state limit. 

Using the above explicit relations between u{t),u{T), 
v^{t) and the nonequilibrium retarded, advanced and 
lesser Green functions, we immediately obtain the tran- 
sient current of Eq. (30) in terms of the nonequilibrium 
Green functions: 

Ia{t) = -^^^l dTTr{s^(i,T)G<(r,i) 

+ S<(t,T)G»(T,t)}. (A.IO) 
This current has exactly the same form obtained form 



the nonequilibrium Green function technique. How- 
ever, as we have pointed out in the practical applications, 
one usually uses the steady-state lesser Green function, 
namely ignoring the first term in Eq. (A.9). This term 
vanishes in the steady-state limit so that it will not af- 
fect the steady-state current. It can also be dropped if 
one takes the initial time to — > — oo so that the cen- 
tral region is assumed to be in an empty state initially. 
However, this term which explicitly depends on the ini- 
tial single particle reduced density matrix (including the 
initial electron occupation in each level and the electron 
quantum coherence between different levels in the cen- 
tral region) is crucial for practical manipulation of a real 
quantum device. Only in the wide band limit (WBL) 
where the non-local time correlation function 



(A.11) 



the integral of the first term in Eq. (A.IO) is reduced to 
the single particle reduced density matrix: G^ {t, t) = 
ip'^^\t)), which results in 



Ia{t) 



dw 

/a(w) 

TT 



X / drIm{e-*l'^(*-^)+^/^''^^'^<'(^')]G"(r,i)}". 

J to 

(A.12) 

Thus, the ignored initial occupation dependence is recov- 
ered in the WBL. 



* Electronic address: wzhangOmail .ncku. edu.tw 

^ D. Goldhaber-Gordon, H. Shtrikman, D. Mahalu, D. 

Abusch- Magder, U. Mcirav, and M. A. Kastner, Nature 

(London) 391, 156 (1998). 
^ S. M. Cronenwott, T. H. Oosterkamp, and L. P. Kouwen- 

hoven. Science 281, 540 (1998). 
^ W. Lu, Z. Ji, L. Pfeiffer, K. W. West, and A. J. Rimberg, 

Nature (London) 423, 422 (2003). 

* T. Hayashi, T. Fujisawa, H. D. Cheong, Y. H. Jeong and 
Y. Hirayama, Phys. Rev. Lett. 91, 226084 (2003). 

^ J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, A. 

Yacoby, M. D. Lukin, C. M. Marcus, M. P. Hanson, A. C. 

Gossard. Science 309, 1280 (2005). 
® T. Fujisawa, T. Hayashi, and S. Sasaki, Rep. Prog. Phys. 

69, 759 (2006). 

R. Hanson, L. P. Kouwenhoven, J. R. Petta, S. Tarucha, 
and L. M. K. Vandersypen, Rev. Mod. Phys. 79, 1217 

(2007). 

* C. Bruder, and H. Schoeller, Phys. Rev. Lett. 72, 1076 
(1994). 

® J. Konig, J. Schmid, H. Schoeller, and G. Schon, Phys. 

Rev. B 54, 16820 (1996). 
^° N. S. Wingroon, A.-P. Jauho, and Y. Meir, Phys. Rev. B 

48, 8487 (1993). 
" A.-P. Jauho, N. S. Wingreen, and Y. Meir, Phys. Rev. B 

50, 5528 (1994). 



Q.-F. Sun, B.-G. Wang, J. Wang, and T.-H. Lin, Phys. 
Rev. B 61, 4754 (2000). 

Z.-G. Zhu, G. Su, Q.-R. Zheng, and B. Jin, Phys. Rev. B 

68,224413 (2003). 

F. B. Anders and A. Schiher, Phys. Rev. Lett. 95, 196801 

(2005). 

X. Q. Li, J. Y. Luo, Y. G. Yang, P. Cui, and Y. J. Yan, 
Phys. Rev. B 71, 205304 (2005). 

M. T. Lee and M. W. Zhang, Phys. Rev. B 74, 085325 

(2006) 

J. Maciejko, J. Wang, and H. Guo, Phys. Rev. B 74, 
085324 (2006). 

J. S. Jin, X. Zheng, and Y. J. Yan, J. Chem. Phys. 128, 
234703 (2008). 

S. Weiss, J. Eckel, M. Thorwart, and R. Egger, Phys. Rev. 
B 77, 195316 (2008). 

T. L. Schmidt, P. Werner, L. Muhlbacher, and A. Komnik, 
Phys. Rev. B 78, 235110 (2008). 

M. W. Y. Tu and W. M. Zhang, Phys. Rev. B 78, 235311 
(2008). 

P. Zedler, G. Schaller, G. Kiesslich, C. Emary, and T. 
Brandos, Phys. Rev. B 80, 045309 (2009). 
X. Zheng, J. S. Jin, S. Welack, M. Luo and Y. J. Yan, J. 
Chem. Phys. 130, 164708 (2009). 

H. Haug and A. P. Jauho, Quantum Kinetics in Transport 
and Optics of Semiconductors, (Springer- Verlag, Berhn, 



17 



1998), springer Scries in Solid-Statc Sciences 123. 

K. Huang, Statistical Mechanics, (2Ed, John Wiley & Sons, 

1987). 

C. Cercignani, Theory and Application of the Boltzmann 
Equation, (Scottish Academic Press, Edinburgh, 1975). 
H. Smith, H. H. Jensen, Transport Phenomena, (Claren- 
don, Oxford, 1989). 

H. J. Carmichacl , An Open Systems Approach to Quan- 
tum Optics, Lecture Notes in Physics, Vol. ml8, (Springer- 
Verlag, Berlin, 1993). 

H.-P. Breuer and F. Petruccione, The Theory of Open 
Quantum Systems, (Oxford University Press, Oxford, 
2002). 

3° J. Schwinger, J. Math. Phys. 2, 407 (1961). 

L. P. Kadanoff, G. Baym, Quantum Statistical Mechanics, 

(Benjamin, New York, 1962) 
^2 L. V. Keldysh, Sov. Phys. JETP, 20, 1018 (1965). 

R. P. Feynman and F. L. Vernon, Ann. Phys. 24, 118 

(1963). 

^'^ K. C. Chou, Z. B. Su, B. L. Hao, and L. Yu, Phys. Rep. 
118, 1 (1985). 

J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 (1986). 
Y. Meir, N. S. Wingreen, and P. A. Lee, Phys. Rev. Lett. 
66, 3048 (1991). 

Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 

(1992). 

^® S. Datta, Electronic Transport in Mesoscopic Systems, 
(Cambridge, England, 1995). 

A. J. Legget, S. Chakravarty, A. T. Dorsey, M. P. A. Fisher, 
A. Garg and W. Zwerger, Rev. Mod. Phys. 59, 1 (1987). 
^° W. H. Zurek, Phys. Today 44 (10), 36 (1991); Rev. Mod. 
Phys. 75, 715 (2003). 

A. O. Caldeira and A. J. Leggett, Pliysica A 121, 587 
(1983). 

F. Haake, and R. Reibold, Phys. Rev. A 32, 2462 (1985). 

B. L. Hu, J. P. Paz, and Y. H. Zhang, Phys. Rev. D 45, 
2843 (1992). 

U. Weiss, Quantum Dissipative Systems, (3rd Ed. World 

Scientific, Singapore, 2008). 

E. Calzetta and B. L. Hu, Nonequilihrium Quantum Field 
Theory, (Cambridge University Press, New York, 2008). 
■^"^ J. H. An and W. M. Zhang, Phys. Rev. A 76, 042127 

(2007) . 

C. H. Chou, T. Yu, and B. L. Hu, Phys. Rev. E 77, 011112 

(2008) . 

J. H. An, M. Feng and W. M. Zhang, Quantum Infom. 
Comput. 9, 0317 (2009). 

T. Yu and J. H. Eberly, Phys. Rev. Lett., 93, 140404 
(2004). 

^° H. Schoeller and G. Schon, Phys. Rev. B 50, 18436 (1994). 
J. Konig, H. Schoeller, and G. Schon, Phys. Rev. Lett. 76, 
1715 (1996). 

M. T. Lee and W. M. Zhang, J. Chem. Phys. 129, 224106 
(2008). 

J. Lehmann, S. Kohler, P. Hanggi, and A. Nitzan, Phys. 
Rev. Lett. 88, 228305 (2002). 
5* J. N. Pedersen and A. Wacker, Phys. Rev. B 72, 195330 



(2005). 

P. Cui, X. Q. Li, J. S. Shao, and Y. J. Yan, Phys. Lett. A 

357, 449 (2006). 
^'^ S. Welack, M. Schreiber, and U. KleinekathLofer, J. Chem. 

Phys. 124, 044712 (2006). 
^"^ U. Harbola, M. Esposito, and S. Mukamel, Phys. Rev. B 

74, 235309 (2006). 

X. Zheng, J. S. Jin and Y. J. Yan, New J. Phys. 10, 093016 

(2008). 

^® L. D. Faddeev, Gauge Fields (Benjamin/Cummings, 1980) 
®° W. M. Zhang, D. H. Feng, and R. Gilmore, Rev. Mod. 

Phys. 62, 867 (1990). 
®^ Usually the stationary path or stationary phase method 
is an approximation to path integral calculations. When 
the path integrals can be reduced to Gaussian integrals, 
the stationary path method will lead to an exact solution. 
For more detailed discussion, see R. P. Feynman and A. R. 
Hibbs, Quantum Mechanics and Path Integrals (McGraw- 
Hill, 1965). 

®^ J. Fransson, O. Eriksson, and I. Sandalov, Phys. Rev. B 

66, 195319 (2002). 
®^ S. A. Gurvitz, Phys. Rev. B 56, 15215 (1997); S. A. 

Gurvitz and Y. S. Prager, Phys. Rev. B 53, 15932 (1996). 

D. Mozyrsky, L. Fedichkin, S. A. Gurvitz, and G. P. 

Berman, Phys. Rev. B 66, 161313 (2002). 
®^ C. Flindt, T. Novotny, and A.-P. Jauho, Phys. Rev. B 70, 

205334 (2004); T. Novotny, A. Donarini, C. Flindt, and 

A.-P. Jauho, Phys. Rev. Lett. 92, 248302 (2004). 

G. Schaller, P. Zedler, and T. Brandes, Phys. Rev. A 79, 

032110 (2009). 

A. Fhijiwara, N. M. Zimmerman, Y. Ono, and Y. Taka- 

hashi, Appl. Phys. Lett. 84, 1323 (2004). 

A. Fujiwara, K. Nisliiguclii, and Y. Ono, Appl. Phys. Lett. 

92, 042102 (2008). 
®^ W. P. Kirk and M. A. Reed (Eds.), Nano structures and 

Mesoscopic Systems (Academic, San Diego, 1992). 
''° M. Switkes, C. M. Marcus, K. Campman, and A. C. Gos- 

sard, Science 283, 1905 (1999). 

M. D. Blumenthal, B. Kaestncr, L. Li, S. Giblin, T. J. B. 
M. Janssen, M. Pepper, D. Anderson, G. Jones, and D. A. 
Ritchie, Nature. Phys. 3, 343 (2007). 
^2 X. Zheng, J. S. Jin, and Y. J. Yan, J. Chem. Phys. 129, 
184112 (2008). 

'■^ S. E. Nigg, R. Lopez, and M. Biittiker, Phys. Rev. Lett. 
97, 206804 (2006). 

J. Wang, B. Wang, and H. Guo, Phys. Rev. B 75, 155336 

(2007) . 

'■^ X. Zheng, F. Wang, C. Y. Yam, Y. Mo, and G. H. Chen, 

Phys. Rev. B 75, 195127 (2007). 
'■^ Y. Mo, X. Zheng, G. H. Chen and Y. J. Yan, J. Phys.: 

Condens. Matter 21, 355301 (2009). 

Y. C. Chang and D. M.-T. Kuo, Phys. Rev. B 77, 245412 

(2008) . 

M. G. Schultz and F. von Oppen, Phys. Rev. B 80, 033302 

(2009) . 

'■^ W. M. Zhang and L. Wilets, Phys. Rev. C 45, 1900 (1992). 



